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Abstract 

We calculate the baryon asymmetry generated at the electroweak phase transition in the min- 
imal supersymmetric standard model, using a new method to compute the CP-violating asym- 
metry in the Higgsino flux reflected into the unbroken phase. The method is based on a Higgs 
insertion expansion. We find that the CP asymmetry at leading order is proportional to the 
change in tan/3 in the bubble wall, which is at most of order 10~ 2 , while at next-to-leading 
order this suppression factor disappears. This result may enhance the final baryon asymmetry 
generated during the electroweak phase transition for small A/? (< 10~ 3 ). 



1 Introduction 
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The baryon to entropy ratio in the observed part of the Universe is constrained to be ns/s ~ 10 
by primordial nucleosynthesis measurements [|I|. Sakharov |2| established more than thirty years ago 
the three basic requirements for obtaining this baryon asymmetry as a result of particle interactions in 
the early universe, namely baryon number violation, C and CP violation and departure from thermal 
equilibrium. These conditions may be satisfied at weak scale temperatures, if the electroweak phase 
transition is first order ||. Electroweak baryogenesis provides an explanation of the observed baryon 
asymmetry of the Universe (BAU) in terms of experimentally accessible physics, hence much attention 
has been devoted to the study of this possibility P|-ffl6|. 

Although the Standard Model contains all the necessary ingredients for electroweak baryogenesis, 
the phase transition is too weakly first order to avoid the wash out of the generated baryon asym- 
metry, for the Higgs mass experimentally allowed [XT]. Moreover, the CP asymmetry induced by the 



Kobayashi-Maskawa phase is far too small to account for the observed Ub/s ratio @, [7|. Therefore, for 
the baryon asymmetry to be generated at the electroweak phase transition, new physics is required at 
the weak scale. 

Among the different extensions of the Standard Model, low energy supersymmetry is a well mo- 
tivated possibility, and thus several groups have recently studied under which conditions electroweak 
baryogenesis is feasible in the framework of the Minimal Supersymmetric Standard Model (MSSM). 
Regarding the strength of the phase transition, a region in the space of supersymmetric parameters has 
been found where the phase transition is strong enough to avoid the wash out of the generated baryon 
asymmetry [I8|-[28|. Such region corresponds to a light Higgs boson and a light top squark, within the 
reach of LEP2 and Tevatron colliders. As for the baryon asymmetry, in the region of parameter space 
favored by the previous phase transition studies, it is mainly generated by charginos and neutralinos, 
provided they are not much heavier than the critical temperature (T c ~ 100 Ge V) |TT|, |13| . 

The physics of the mechanism which produces the baryon asymmetry is essentially agreed upon, 
namely particles in the plasma interact with the bubble wall and due to CP-violation lead to a chiral 
asymmetry of fermions in front of the wall, which in turn biases sphaleron processes to produce 
the BAU. Several groups have estimated the baryon asymmetry generated at the electroweak phase 
transition in the MSSM: in refs.flll], CP-violating source terms were computed and then inserted 
into a set of diffusion equations, with a prescription not motivated by first principles. Subsequently, 
the correct definition of the CP- violating sources was self-consistently derived in [16j , using the closed 
time-path formalism to write down a set of quantum Boltzmann equations describing the local particle 
densities. The CP-violating sources were computed using a mass expansion |TT[ and a Higgs insertion 
In 



expansion 
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the diffusion equations and the source terms were derived together within 
the WKB quasi-particle approximation. 

TT| , |13| , m, the source terms are proportional to V2{z)d z v\(z) — v\(z)d z V2(z), which is zero for 

8f that the angle j3 varies at most 



In 

tan (3 



v%jv\ constant in the bubble wall. It has been shown |29], 
by a few percent over the wall in the MSSM, thus it results in an important suppression factor. On 
the other hand, in [TT| it was mentioned that for charginos and neutralinos such A/3 dependence could 
disappear at higher orders in the mass expansion used to calculate the source terms, and in [TJj], CP 
asymmetries of the reflection and transmission amplitudes for charginos are computed numerically, 
and they are non-zero for constant tan (3 []. 

1 The authors of pjJ] also find a non- vanishing Higgsino source term for constant tan/3, but they consider a different 
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Motivated by this discrepancy and the possibility that it is a consequence of the lowest order 
approximation, we have performed a new computation of the baryon asymmetry generated by charginos 
at the electroweak phase transition within the MSSM. Although we introduce a different method to 
compute the CP-violating source terms, it is based on a Higgs insertion expansion similar to the one 



employed in [T3|, [LJ|, thus we expect that our conclusion will also apply to that approach. 

The remainder of the paper is structured as follows. In section 2 we describe the method used to 
compute the CP-violating asymmetries, and present the results for Higgsinos at leading and next-to- 
leading order in the Higgs insertion expansion. In section 3 we compute the baryon number induced 
by the Higgsino current, and we conclude in section 4. 



2 CP asymmetry for Higgsinos 

Baryogenesis is fueled by CP asymmetries induced by the advancing bubble wall. Unremovable CP-odd 
phases appear in the mass matrices due to either 

a) CP-violating interactions in the thermal loops that correct the dispersion relations of particles 
propagating in the plasma [Q. 

b) space-time dependence of the scalar vevs inside the bubble wall (for more than one Higgs field), 
which induces space-dependent CP-violating phases. These phases cannot be rotated away in two 
adjacent points by the same unitary transformation 



When present, the second mechanism dominates over the first one, since in the first mechanism 
there are suppression factors coming from loops. While in the Standard Model only the first mechanism 
is possible, the second one controls the generation of the baryon asymmetry in all the extensions of the 
Standard Model proposed in the literature for electroweak baryogenesis. Thus, particle mass matrices 
acquire a non-trivial space dependence when bubbles of the broken phase nucleate and expand during a 
first order electroweak phase transition. This provides fast non-equilibrium CP-violating effects inside 
the bubble walls (thick walls) or in front of the bubble walls (thin walls), and may generate a baryon 
asymmetry through the anomalous (B + L) -violating processes when particles diffuse to the exterior 
of the bubble, in the unbroken phase. 

In what follows, we will focus on CP-violating effects in the Higgsino current, since they make 
the dominant contribution in the region of the MSSM parameter space preferred by phase transition 
studies jn|. However the method presented here is completely general, and may be applied to any 



particle whose mass depends on the scalar vevs. We will perform the calculation in the thin wall 
regime (L w < £, where L w is the bubble wall width and I the mean free path of the particle). This is a 
reasonable approximation for weakly interacting particles, such as charginos. For bubble walls thinner 
than the mean free path, the incoming fermions interact with the bubble wall like quantum mechanical 
particles scattering from a potential barrier. CP-violating interactions with the scalar field result in 
different reflection probability for fermions of a given chirality and their corresponding antifermions, 
leading to a CP asymmetry in the reflected chiral number flux [Q. 

The chargino mass matrix is, in the basis of Winos and Higgsinos (W, H) , 




m 

v gvi/V2 

Higgsino density |3(J , namely the difference of the two helicity states instead of the sum, which is the one considered by 
all the other groups, including us. 
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with Vi the spatially varying Higgs field vevs and m, fj, the soft supersymmetry breaking parameters, 
which contain one physical CP violating phase (/i = |/z|e^). In the symmetric phase (t>, = 0), Winos 
and Higgsinos are the mass eigenstates, while in the bubble wall and the broken phase they mix to 
form the mass eigenstates, denoted x% > i = 1? 2. We treat the wall as planar, and assume it has reached 
a static configuration in the wall rest frame, with the vevs of the scalar fields being functions only of 
the z coordinate. 

The quantity of interest is the CP-violating asymmetry in the Higgsino current reflected into the 
unbroken phase. In the rest frame of the bubble wall, it is given by 

2 

jcp = iU + iLf + Ei^ ( 2 ) 

i=i 



+ / ^3 EO^fe.rf - \T^ 6 (pi,P,)\ 2 ) ^ P \E, Pl ), (3) 

where 

p m = ^E 2 -p\-m 2 , p^ = sJe? ~p\- /i 2 , Pi = \JE 2 ~ P\\— m l ( 4 ) 
are the absolute value of the z component of the Wino, Higgsino and chargino momenta, and 

P u (E, Pz ) = e(E+v J z)/T + 1 , p b (E,p z ) = e(E _ VwPz)/T + 1 , (5) 

are the thermal distributions of the charginos in the unbroken (u) and broken (b) phases, as seen from 
the rest frame of the wall. v w is the wall velocity, v w ~ 0.1 in the MSSM. The use of equilibrium particle 
distributions is a good approximation, because any departure from thermal equilibrium is caused by 
the passage of the wall, and therefore is 0(v w ). Since we will see that the final CP asymmetry is 
already linear in v w , working with thermal equilibrium distribution functions amounts to ignore terms 
of higher order in v w [II]. 

It is well known that one-loop self-energy corrections to the propagator modify the dispersion 
relations of the particles in the plasma. The main thermal effects are that particles propagating in 
the plasma acquire an effective mass (even if they are massless in vacuum) and have a finite life-time 
(damping), due to incoherent thermal scattering with the medium. The damping rate, 7, is defined as 
(minus) the imaginary part of the solution uo = uj(k) of the dispersion relation. 

For the Higgsino and Wino, the effective plasma masses in the thermal bath may be well approx- 
imated by their value in the present vacuum (m^-(T) ~ |/i| 2 ^^(T) ~ m 2 ). However, it has been 
shown that the effects of damping can lead to a sizeable suppression of the CP asymmetry || 0, due 
to the loss of coherence of the wave function. The damping rate of Winos and Higgsinos has been 
estimated in |H| to be 7^ ~ 0.025T ,j w ~ 0.065T, hence the mean free path £ ~ 1/(27) ~ (10-20)/T 



is comparable to the wall width L w ~ (20 — 30) JT f2"9| , and decoherence effects may be relevant. 

The effects of damping may be taken into account by including the imaginary part of the fermion 
self-energy in the dispersion relation || 0, |13|]. In our case, this leads to the approximate dispersion 
relation 

[cu(k)+2 7 ] 2 = k 2 + m 2 . (6) 
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If we choose u to be real, the momenta must become complex in order to satisfy the dispersion 
relations, and propagation of particles in space is damped. We have taken u to be real because 
energy is conserved in the scattering off the wall, so the reflection and transmission probabilities are 
time independent. We could have satisfied the dispersion relations with real momenta and complex 
lu, but then the reflection and transmission amplitudes would have an exponentially decaying time 
dependence, which would require us to study the time and space dependence of the particle scattering 
process. To obtain an estimate of the damping effects, from (§) we approximate 

k ~ ± < \JuS 2 — m 2 + i'-f i ~ ±(V uj 2 — m 2 + ry) , (7) 

( yuj 2 — m 2 J 

where k = |k|. 

We thus calculate the reflection and transmission amplitudes at zero temperature (except for the 
damping), using the LSZ reduction formulae in terms of the propagator in the presence of the bubble 
wall H: 



A = J d 4 xj d 4 y e' iqiX e iq f y u{q f ){id - fi)S(y, x)(-id - m)u(ft) 

= (27r) 3 %? - qf ) S(qj - qf) 8(E f - Eft A(qf, qf), (8) 

with 

%,x) = (0|T[*(y)*(z)]|0). (9) 

An analogous expression holds for antiparticles. The spinors in formula @ are on-shell and normalized 
to unit flux in the z direction, i.e. 

u j z u = 1. (10) 

Momenta in the x and y directions are conserved, because the potential created by the bubble 
wall only depends on the z coordinate. The transmission and reflection amplitudes are then functions 
only of the momenta in the z direction and can be computed in a simpler way by first boosting to a 
frame where q x ,q y = 0. With the proper normalization chosen for the spinors (p!U|), the amplitude in 
the boosted frame is simply given by with the propagator and incoming and outcoming momenta 
substituted by the boosted ones. 

The expression of jcp can be further simplified by using CPT and unitarity constraints, which 
imply 

i^i 2 + i^-^r + e \T b ^\ 2 = 1 . (ii) 

i=l 

Substituting (|TT|) in eq.@, and expanding the Fermi distributions for small wall velocities we obtain: 



jcp = yJ-^e {Hp»+pi)[\R^ 6 (-p»,pX-\R u h- 



|2l 



(Pm + Pi ) [ I Rjv^ff ( ~Pm , Py, ) 1 2 - I Rw->H ( ~Pm , P» ) 1 2 1 



+ (p 2 - Pl )[|T^(p 2 ,^)| 2 -|f^(p 2 ,^)| 2 ]}^p(E) [l-p(E)] , (12) 

where p(E) = \j{e E l T + 1) is the Fermi distribution. This result explicitly shows that the out of 
equilibrium condition needed for baryogenesis is due to the expansion of the bubble wall through the 
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thermal bath. As anticipated, the CP violating current jcp is linear in the wall velocity, v w , and the 
use of the equilibrium distributions in the calculation is justified. 

In order to compute the amplitudes in eq.(|8|) we would need the exact propagator in the presence 
of the wall (H). In ref. [11]], the approach was to perform an expansion in powers of mass, which 



is effectively an expansion in M(z)/E, but this approximation is not justified, since the region of 
interest is always E ~ M. Instead, we perform a Higgs insertion expansion |13|], which will be a good 
approximation at least close to the symmetric phase. By making a phase redefinition of the Higgsino 
field, we can write the mass matrix as 

M x (z) = M° x + 5M x (z) (13) 

where 

M ° = ( o M ) (14 » 

and 

with Ui(z) = gvi(z) / \/2. Now, expanding in SM x (z), the approximate result for the propagator is 

S(x 2 , x x ) = [ J] dzi S^(x 2 , zx)6M x (zx)S^(zx, z 2 )6M x (z 2 ) . . . S^(z u Xl ), (16) 

where the integration is done over all Z{(— oo, oo) and stands for the propagators of Higgsino and 
Wino in the symmetric phase, at zero temperature. When the approximate dispersion relation (|7|) is 
used, includes also the damping. In the boosted frame (p x — p y — 0), for the Wino it is given by 

c(0), \ . fdp z e iP*(z2-*i) f . 

S {0 \z 2 , Zl ) = — (u 7 °-p gl * + Tr A 

2 { \Pm Pm) 

+ 9(zx - Z2 ) e -(^-7#)(^-«) (^ 7 o + ^+-U . (17) 

\Pm Pm J ) 

The Higgsino propagator has an analogous expression, changing m — > p m — > p^ and 7^ — > 7^. 



Here, u is the energy in the boosted frame, u = \jp m + m? = ^jp 2 ^ + |/i| 2 and we have only included 
the damping rate in the exponential factors, where the effect is expected to be more important. 

2.1 Leading order 

The leading CP violating contributions arise at second order in the Higgs insertion expansion, 0(v 2 ), 
as expected from refs. |TT1, |H| where this was also the order of the CP violating currents. Since the 
perturbation matrix 8M x (z) is off-diagonal in the Wino-Higgsino basis, at lowest order only Winos 
can be reflected or transmitted into Higgsinos, and we get 

R$^ 6 (-Pm,Pv) = %(?/) r dze^ + ^ z e-^{ Ul {z)L + u 2 (z) e^R}u w ( qi ) , (18) 

J 
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where 7 = 7^ + 7^, q { = (u,-p m - 27^), g/ = (u,Pp + 27.fr)> oj is the energy in the boosted 
frame (q x = q y = 0) and L, R are the chiral projectors. The result for the transmission amplitude 

^Lh(p™>ZV) can De obtained just changing q { by g. = (w,p m - £7^) and p m by -p m in eq.flTJ). 

Using standard techniques to compute Dirac matrix traces, and substituting the leading order 
reflection and transmission contributions in eq. (|T2"D , we obtain 

jgl = ^ sin0| F <yl \p m ) p(E) [1 - p{E)\ , (19) 

where F^(p m ) = -(p m + p^)FP(p m ) + (p m -p»)Ff\p m ), being 

Fi x \p m ) = / dz 1 dz 2 e l{p - +p ^- Z2) e-^ zl ~ Z2) {u 1 (z 1 )u 2 (z 2 ) - u 1 (z 2 )u 2 (z 1 )} (20) 

the contribution coming from the reflection, and 

F?\p m )=F?\-p m ) . (21) 

the transmission one. ^From the above equations, we see that the CP asymmetry j^p vanishes if 
tan/3 = u 2 /u\ is constant along the bubble wall, in agreement with refs. JTT , |T^| , where the computation 
of the CP violating currents was done to lowest nontrivial order in the mass and Higgs insertion 
expansions, respectively. 

To obtain the final result we must still specify the shape and speed of the bubble wall. For the 
wall profile we take the following semirealistic approximation | : 

v(z) = ^l- CO s^}i6(z)-6(z-L w )}+v6(z-L w ) 

P(z) = ^ {1 - cob [6{z) -e(z-L w )]+ A/3 Oiz-Lv) (22) 

where v 2 = v \ + v\ and are the vevs of the two Higgs doublets in the broken phase, at the critical 
temperature. For the remaining parameters we take the bubble wall width L w = 25/T, the wall 
velocity v w = 0.1 and Af3 = 0.01 [g, |28| . 



Fig. 1 shows Fj; 1 \p m ), which is the dominant contribution, as a function of p m: for v/T c = 1 and 
different values of the Higgsino and Wino masses, (|/i|,m). The weak coupling at the phase transition 
temperature is a w = 0.035. We can see an enhancement of the CP asymmetry for values of |/i| close 
to m, i.e., when the Higgsino and Wino are nearly degenerate. This resonant behavior, also found in 



13] , is enhanced by the approximations made in our calculation, because the masses of the particles 
propagating in the plasma are those of the symmetric phase, and thus strong degeneracies occur when 
\fi\ ~ m. 

Regarding the damping effects, as explained before the lifetime of the charginos £ ~ 7 _1 is of the 
same order of magnitude as the wall width, L w . In this situation, a priori it is not clear whether 
the damping will have an important effect or not. To estimate it, we have repeated the leading order 
calculation neglecting the damping rate, i.e. setting 7 = 0. The resulting j C p is typically one order of 



2 The difference in the final result using the ansatz (^2|) or the functional kinks is typically 5%, although for some 
values of the soft masses may be as large as 40%. 
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Figure 1: (p m ) as a function of p m for Higgsino mass = 1. and Wino mass m = 0.5 (solid), 
m = 1.5 (dashed-dotted), m = 2. (dashed) and m — 3. (dotted). All in units of the temperature. 

magnitude larger, independently of the Higgsino and Wino masses (provided they are not degenerate 3 ). 
Therefore we conclude that although the suppression of the CP asymmetry due to the scattering of 
charginos with the plasma is not as enormous as in the case of strongly interacting particles (such as 
quarks in the Standard Model), decoherence effects are not negligible. 

2.2 Next-to- leading order 

In this section, we will compute the CP violating current jcp at next order in the Higgs insertion 
expansion. Now both Wino and Higgsino can be reflected or transmitted into Higgsino and the CP 
asymmetry ([12]) can be written as 



where the function F^(p m ) is given in Appendix A. 

Contrary to the lowest order result, we find that y CP ^ even if tan/3 remains constant in the 
bubble wall. Thus the A/3 suppression disappears at higher order in the Higgs insertion expansion, as 
mentioned in [|ll|], and in agreement with the numerical computation of [|l4fl . 

For constant v 2 /vi we obtain F (2 \p m ) = -(p m +p tl )F^ 2) (p m ) + (p m - p^Ff' \p m ), with 



3 For \/j,\ ~ m, since the transmission contribution behaves as l/(p m — Pp) the computation without damping rate is 
divergent. 




(23) 
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Figure 2: Fj; 2 \p m ) as a function of p m for tan (3 = 2, |/i| = 1. and m = 0.5 (solid), m = 1.5 (dashed- 
dotted), m = 2. (dashed) and m = 3. (dotted). All masses in units of the temperature. 
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(24) 



and 



sin(4/3) Im( r dz^z^e^e-^ 4 
2p )l Uo 

dz 3 u(z 3 )e ipsZ:i r dz 2 u(z 2 )e- iprZ2 H dz^z^e'^ 1 
o Jo Jo 



(25) 



where u(z) = gv(z)/ y/2, p s = p m + p M , p r = p m - p^ and 7 = 7^ + j w . 

Again, the dominant contribution comes from the transmission amplitude. In Fig. 2 we plot 
(p m ) as a function of p m , for the same values and shape of the wall profile as in the previous 
section, except that now A/3 = and we take tan/3 = 2 (favored by studies of the phase transition). 
Recall that for this choice of parameters, the CP asymmetry at lowest order vanishes, so this is the 
leading order result. We also find an enhancement of the CP asymmetry when |/i| ~ m. 



3 Baryon asymmetry 



Now we have to solve a set of coupled differential equations describing the effects of diffusion [32 



particle number changing reactions and CP asymmetries, to obtain the various particle densities in 
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the MSSM []. We follow the approach of refs.JTT], [L3[] (see [[Tl|] for details). If the system is near 
equilibrium and the particles interact weakly, the particle number densities are given by 

T 2 

rii = hHi— , (26) 

where /ij is the local chemical potential for particle species i, and hi is a statistical factor of order 2 
(1) for light bosons (fermions), while for particles much heavier that T it is Boltzmann suppressed. 

The particle densities we need to include are the left-handed doublet, q^ = (t L + b L ), the right 
handed top quark t = tn, the Higgs particles h = (H® + + + H^), and their superpartners, 
QL,tR, h. The particle numbers of these species can change due to the top quark Yukawa interaction 
with rate T t , the top quark mass interaction with rate T m , the Higgs self-interactions with rate Th, 
the strong sphaleron interactions with rate T ss , the anomalous weak interactions with rate r„, s and 
the gauge interactions (which we shall assume that are in equilibrium). Then, the system may be 
described by the densities Q = q + q, T = t + t and H = h + h. Alternatively, one could use the Higgs 
particle density H' = (H® + H{ + H$ + H%)+ superpartners, which does not have the A/3 suppression 



in the corresponding source term (T^|. However, if Yukawa or helicity-flipping interactions were in 
equilibrium this would force Hi = —H 2) leading to H' — Q = T = 0. Thus, the approximations 
we make below would give a vanishing baryon asymmetry, and a numerical solution of the diffusion 
equations for the three particle densities including all the relevant interaction rates is needed j30 |. 



CP-violating interactions with the phase boundary produce an injected Higgsino flux, which we 
model as 

J inj (z)=Ucp6(z-v w t) , (27) 

where jcp is the net Higgsino flux reflected into the unbroken phase and £ defines the persistence 
length of the current in the vicinity of the wall, i.e., it parameterizes our ignorance about how the 
injected flux thermalizes. This approximation is reasonable if the injected current thermalizes in a 
time r t h short in comparison to the time the particle spends diffusing before being recaptured by the 
wall, i.e., for small velocities of the wall (v w l/v^) [P^H - We use the estimate for £ of ref.jjlO 



£ ~ 6Dh(v), where is the diffusion constant of the Higgsino, which can be approximated by the 
one of left handed leptons, ~ 110/T fllU| , and (v) is the average velocity of the Higgsinos in the 
reflected flux, 

^ J$*F®(Pm)%p(E) [I-P(E)] 
f^F(H Pm )±p(E)[l-p(E)} ' 

with F®(p m ) (i = 1,2) defined in eqs. (|T9|) and (p3|), respectively, for the leading and next-to-leading 
order computations of the CP asymmetry. 

In ref. the closed time-path formalism was used to derive a set of quantum Boltzmann equa- 
tions describing the local number density asymmetries of the particles involved in supersymmetric 
electroweak baryogenesis. In these diffusion equations the CP-violating sources which fuel baryogene- 
sis are self-consistently incorporated. According to them, the CP-violating source term which should 
be inserted in the diffusion equation for the Higgs density is given by 

jinj 

T h , (29) 

T 



4 We neglect the Debye screening of induced gauge charges, since the effect on the baryon number produced is 0(1) 
53l. 



9 



where r = 7^. is the thermalization time of the Higgsino. 

Assuming that the rates T t and T ss are fast, so that Q/kQ — H/k H — T/k T = 0(1/T t ) and 2Q/kg — 
T/k T + 9(Q + T)/k B = 0(l/r ss ), we obtain 

« = ^f!V!t , +0(i/r„,i/r.), 

k H (9k T + k B + 9k Q ) 

T = -» , wfi\ +9 ft i + °(Vr,, i/r.) • (30) 

k H {9k T + k B + 9 k Q ) 

Substituting these expressions we find the equation for the Higgs density 

DH" - v w H' - TH + 7 = , (31) 

where D is an effective diffusion constant, T is an effective decay constant and 7 is an effective source 
term, given by Q 



D 



D q (9k Q k T + k Q k B + 4k T k B ) + D h (9k T + k B + 9k Q )k H 
9k Q k T + k Q k B + 4k T k B + k H {9k T + k B + 2fc Q ) 



r = rr + r 1 9fc T + fc B + 9fc Q f . 

1/1 mJ 9k Q k T + k Q k B + 4k T k B + k H (9k T + k B + 2k Q ) ' 1 ; 

/c H (9/c T + k B + 2k Q ) 

1 = T h 



9k Q k T + k Q k B + Ak T k B + k H (9k T + k B + 2/c Q ) 

with D q (Dh) the diffusion constant for quarks and squarks (Higgs and Higgsinos). 

We are interested in an analytic solution to eq.(|3T|) which satisfies the boundary conditions 
if (±00) = 0, and at the interphase z = 

H\± = 0, DH'\ + - = --]cp , (33) 

T 



which are derived by integrating up the eq. (plf) through z — 0, imposing the condition that H is 
at most step-like discontinuous across the wall. Here, jcp = (7/7^) jcp- We use a z— independent 
effective diffusion constant and a step function for the effective decay rate T = T6(z). The values of 
D and f depend on the supersymmetric parameters, for the considered range D ~ 0.8GeV _1 , f ~ 1.7 
GeV H . 

Then, the solution of eq . ([31]) in the symmetric phase (z < 0) reads 



H(z) = Ae ZVw/D , (34) 

with 



A = ^ (35) 



and A ± = !^g±gg 



3 Our expressions slightly differ from those of ref.[[TT], but we have checked that this difference is numerically negligible. 
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In ref.]Tl|], average current densities J(z) were computed in each point z of the bubble wall, and 



CP- violating sources 7q(z) associated with these currents constructed as 

lQ (z) ~ J\z)/r . (36) 

In order to qualitatively understand the dependence of the produced baryon asymmetry on the various 
parameters, the sources were approximated as step functions of width L w . In the thin wall limit 
(L w — > ,L W J° = constant) the expression of the coefficient A obtained in coincides with ours 



with the substitution £jcp —> L W J°. This result may help to understand the relation between the two 
approaches, in the regime when both are applicable. 

^From the form of ( |3~5]) we see that the CP- violating densities are non zero for a time t ~ D/v^, 
thus the assumptions which lead to the equation ([31]) for the Higgs density are valid provided T t , T ss 
vl/D. 

What we aim to compute is the total baryon number density left inside the bubble, so now we 
need to consider the effect of the weak sphaleron processes, since they provide the only source for net 
baryon number. The equation satisfied by the baryon number density, ub is 

D q n" B - v w n' B - 6(-z)n f T ws n L = , (37) 

where is the total number density of left handed weak doublet fermions, n/ = 3 is the number of 
families and for the weak sphaleron rate we take T ws = Qua^T (k ~ 1). It has recently been estimated 
as T ws ~ Ca w log(l/g)(a w T) 4 f34|, but lattice measurements of the rate are consistent with C ~ l/a w 



35] so this does not affect the numerical value of our result. 

Assuming that all squarks except t^ R and bi are heavy, there is no suppression due to the strong 
sphaleron and n B is given by fll|| 

9kqkT — 8kxkB — 5kqk B 27 
UL= k H (9k T + k B + 9k Q ) H = 82 H - (38) 



Substituting the Higgs density (|3|) in eq. (|37D , we obtain 

n B siADr* 



(39) 



where s = 2tt 2 g*T 3 / '45 is the entropy density, with ~ 126 the effective number of relativistic degrees 
of freedom. 

In order to compare our approach with previous related work, we have calculated the final baryon 
asymmetry ( |3"9"D using both, the leading and the next-to-leading order results for the CP asymmetry 
in the Higgsino current. When we use the lowest order CP asymmetry we obtain a baryon number 
to entropy ratio of the same order of magnitude as in 0, where the CP-violating sources were 



self-consistently incorporated in the diffusion equations. In previous calculations [11, [L3|, different 



definitions of the CP-violating currents were used without any first principles justification, leading to 



smaller numerical results for the final baryon asymmetry (about one order of magnitude in JTI and 
two orders of magnitude in |13|). 

Fig. 3 shows our leading order results. ^From eq.(^), we see that the baryon number produced 
depends linearly on the phase of the Higgsino mass parameter, 0. So in Fig. 3 we plot the value 

6 We thank A. Riotto for pointing out this fact to us. 
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Figure 3: Contour plot of |sin</>| in the plane (|/i|,m) for fixed n^js ~ 4 x 10 -11 , according to our 
leading order result for the CP asymmetry in the Higgsino current. The masses are in units of the 
temperature. 



of | sin 0| necessary to obtain %/s ~ 4 x 10~ n , in the \\x\ — m plane, using the leading order CP- 
asymmetry computation with A/3 = 0.01. The shape of the curves is similar to the ones obtained 
in [13] but we need a smaller phase to explain the observed baryon asymmetry, as discussed above. 
We have also calculated the baryon number to entropy ratio at leading order for other values of the 
parameter A/3, namely A/3 = 0.001, 0.005, an we find that it scales as n B /s oc A/5. 

For tan (3 constant in the bubble wall, the CP asymmetry at lowest order vanishes, and at next-to- 
leading order we obtain 

— ~ (lO- n -lO- 8 )sin0 , (40) 

s 

where the largest values correspond to nearly degenerate Wino and Higgsino. In most cases, | sin0 | 
should be of order 0.1-1 to explain the observed baryon asymmetry. These results correspond to 
tan/3 = 2, but they can easily be converted to other values of tan/3, since the (3 dependence is just 
a global factor (see eqs.(^) and ([25"])). The values above seem to indicate that next-to- leading order 
corrections are relevant only for very small A/3(< 0.001). 



4 Conclusions 

We have computed the baryon asymmetry generated at the electroweak phase transition in the MSSM. 
The leading CP asymmetry which fuels baryogenesis is in the Higgsino current, provided that Higgsinos 
and gauginos are not much heavier than the electroweak critical temperature (T c ~ 100 GeV). The 
main motivation for the calculation was to settle the question of whether or not the A/3 dependence of 
the CP-violating Higgsino current was a consequence of the approximations used in the computations 
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IT], P|, ff§|. Given that A/3 is at most 2.5 x 1(T 2 in the MSSM [gR E3, such dependence leads to a 



significant suppression of the produced baryon asymmetry. 

We have done the calculation in the thin wall regime, using a new method to compute the CP 
asymmetry by expanding in the Higgs vevs. At leading order, our results show the same parametric 
dependence as previous related approaches [II , |T3|| , in particular the CP asymmetry vanishes if tan (3 = 
Vijv\ is constant along the bubble wall. However, at next-to-leading order we find a non-zero CP- 
asymmetry even if tan/3 is constant, in agreement with the numerical calculation of ||14j| . Although 



our approach differs from those of refs.||iT|, [U|] , all of them contain a Higgs insertion expansion and 
thus we believe that this result will also apply to their methods. 

We have estimated the effects of damping in the CP asymmetry, using a simple model for deco- 
herence. In the leading order calculation we find a suppression of order 0.1 due to the incoherent 
scattering of charginos with the particles in the plasma. 

The subsequent solution of the diffusion equations using the leading order CP asymmetry leads to 
a baryon number to entropy ratio in agreement with observation, even for small CP-violating phase 
of the [i parameter, | sin0| ~ (10~ 2 — 10 -4 ). These results are comparable with those of ref.[16[], for 
the same range of the relevant MSSM parameters. At this order in the Higgs insertion expansion, the 
dependence of n^js with A/3 is approximately linear. 

For tan (3 constant in the bubble wall, the CP asymmetry at lowest order vanishes, and the next- 
to-leading order calculation gives a baryon asymmetry sufficiently large if the CP- violating phase is in 
the range | sin0| ~ (0.1 — 1). Smaller values are acceptable only in the region m ~ |/i|. 

The various approximations made in the current analysis lead to a sizeable uncertainty in the final 
result. Since the mean free path for charginos I ~ L w (bubble wall width), we expect the thin wall 
limit to give an 0(1) estimate of the true solution. Regarding the Higgs insertion expansion, our next- 
to-leading order calculation shows that the leading order result is reliable, except for very small values 
of A(3(< 0.001). We have also made a number of simplifications of the diffusion equations [O, which 
are valid provided T t , T ss 3> v^/D, T ws <C f 2 / D and the scattering processes due to Yukawa couplings 
other than top are slow. Given our poor knowledge of the parameters involved (diffusion constants 
D, reaction rates T ttSS ,ws and wall velocity v w ) it is not possible to ensure that these inequalities are 
satisfied, but the present estimates seem to indicate so. 

Our most severe approximation was the insertion of the CP-violating source in the diffusion equa- 
tion for the Higgs density (|31]). We have modeled the injected Higgsino flux as a delta function localized 
at z = v w t. The largest uncertainty comes from our ignorance about how this injected flux thermalizes. 
Following [ 10 1 , we have parameterized it by £, the persistence length of the injected current, estimated 
to be £ ~ 6Dh(v), where GD^ is the velocity randomization time of a diffusing Higgsino and (v) is the 
average velocity of the injected flux. To calculate this parameter more precisely involves going beyond 
the diffusion approximation to determine exactly how the reflected asymmetry enters in the diffusion 
equation fll6|. Taking into account the above uncertainties and approximations, we estimate that our 
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results for the baryon asymmetry are reliable to about one order of magnitude. 

Finally, we want to comment on the differences between our approach and the one in refs. |I3" 
since both are based on a Higgs insertion expansion. Our approach is only applicable in the thin 
wall regime, and we perform a tree level computation of the reflection and transmission amplitudes to 

On the other hand, the method of 



obtain the Higgsino current, much as in [11 



13, O] is valid for any 



thickness of the bubble wall and they calculate directly CP-violating currents using the closed-time 
path formalism, which involves a one-loop computation at finite temperature. However, at least in 
the thin wall regime, both approaches should be somehow related, and it will be very interesting to 
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understand the connection between them. 
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Appendix A. 

We present here the different contributions to the function F^ 2 \p m ) which appears in the CP asym- 
metry at next-to-leading order in the Higgs insertion expansion, (^). Taking into account 
the unitarity relation (|Tl~D, we only need the contributions from the reflection and transmission of 
Winos into Higgsinos, and the reflection of Higgsinos into themselves. We work in the boosted frame 
— % = 0). In the reflection, the momenta of the initial particles are = (u,—p^ — ijjj) for 
Higgsino and % = (u>, —p m — ilw) f° r Wino, while qf = (u;,p M + z7#) is the momentum of the Higgsino 
reflected into the symmetric phase. 
[From H — > H reflection we obtain: 



F£{p m ) = -2p„F^(p m ) 

= —2im\fi\ J dzi . . . dz^ J 



dp z dp' e^ PM_7 ff^ 2l+22 ^e^~ ip ''~ 7 ff^ Z3+Z4 ^e _ * Pz ^ 2 ~^ 1 ^ +ip ' z ^ 4_Z3 ' > 
2n 2n (p 2 — m 2 + i"iy/){p' 2 — m 2 — ijp) 

x {2m\p,\ cosi^^Owi^W^W^) - u 1 (z 1 )u2{z2)u 1 (z 3 )u 2 {z4)) 

+(q f -p )[U2(ZI)U 1 (Z2)U2(Z 3 )U 2 (Z4) ~ Ui(z 1 )u 2 (z 2 )u 1 (z 3 )u 1 (z4}} 
+ (Qi ■ P)[U2(Z 1 )U 2 (Z2)U 2 (Z 3 )U 1 (Z 4 ) - Ui (Zi)Ui(z 2 )Wi (23)^2(24)] 
+ (q f ■ p)[ui{zx)ui{z2)u 2 {Zz)ui(Zi) - U 2 {zi)u 2 {z2)u 1 {z 3 )u2{z 4 )] 

+ {Qi ■ PlMz^Ui^Ui^Uiiz^) - U 1 (Z 1 )U2{Z2)U2{Z 3 )U 2 (Z 4 ))} , (A. 1) 

where p 2 = uo 2 — p 2 ,, p' 2 = u 2 — p' 2 . By using symmetry arguments, it is easy to see that the above 
expression vanishes identically. 

The contribution coming from W — *> H transitions can be written as 



F w(Pm) = "(Pm + P^tiPm) + (Pm ~ P,)F^%( Pm ) , (A. 2) 



where 
F (2) - (v ) 

1 rW^ m > 



2m\n\ ( f f dp z dp' z e^ Ptl 7 » 



Im y J dz\ . . . dz4 J 



(*PM-7&)23 + («Pm-7w) z lp[- i (PM+Pm)-7]^4piP2(2;3-Z2)+ip / z (22-2l) 



Pv V J J 2tt 2tt (p 2 - m 2 + ^7iy)(P' 2 ~ + ^Ih) 

x {2m|/i| cos(p[u2(z 1 )u 1 (z 2 )u2(z 3 )u 1 (z 4: ) - u 1 (z 1 )u 2 (z2)u 1 (z 3 )u 2 (z4)} 
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+ (<?/ • qi)[u2(z 1 )u 1 (z 2 )u2(z 3 )u 2 (z 4 ) - U 1 (Z 1 )U 2 (Z 2 )U 1 (Z 3 )U 1 (Z 4 )] 
+ (<?/ ■ P)[u2{zi)u 1 {z 2 )u 1 {z z )Ui{z i ) - U 1 (Z 1 )U2{Z 2 )U2{Z 3 )U2{Z4)] 
J r{qi-V')[Ux{Zi)u 1 {z 2 )u2{z^)u 1 {Zi) - U 2 {Z 1 )U2(Z2)U 1 (Z 3 )U 2 {Z II )) 

+ (p ■ P')[u2(zi)u2(z2)u 2 (z 3 )u 1 (z4) - U\ (z x )u\ {z 2 )ux {z 3 )u 2 (^4)] }) , (A. 3) 

with 7 = 7^ + The contribution of the transmission from the broken phase, F^(p m ), can be 
obtained from eg. ( [A. 3|) changing — > q\ = (uj,p m — i"fw) and p m — > — p m in the exponentials. We 
have only included the damping rate in the exponential factors. 

References 

[1] G. Steigman, Ann. Rev. Astron. Astrophys. 14 (1976) 339 
[2] A.D. Sakharov, JETP Lett. 6 (1967) 24 

[3] V.A. Kuzmin, V.A. Rubakov and M.E. Shaposhnikov, Phys. Lett. B155 (1985) 36 

[4] For a review, see A.G. Cohen, D.B. Kaplan and A.E. Nelson, Ann. Rev. Nucl. Part. Science 43 
(1993) and references therein. 

[5] G.R. Farrar and M.E. Shaposhnikov, Phys. Rev. Lett. 70 (1993) 2833; Phys. Rev. D50 (1994) 
774 

[6] M.B. Gavela, P. Hernandez, J. Orloff, O. Pene anc C. Quimbay, Mod. Phys. Lett. A9 (1994) 795; 
Nucl. Phys. B430 (1994) 382 

[7] P. Huet and E. Sather, Phys. Rev. D51 (1995) 379 

[8] D. Comelli, M. Pietroni and A. Riotto, Phys. Lett. B354 (1995) 91; Phys. Rev. D53 (1996) 4668 

[9] M. Joyce, T. Prokopec and N. Turok, Phys. Lett. B338 (1994) 269; Phys. Rev. Lett. 75 (1995) 
1695; erratum, ibid. 75 (1995) 3375; Phys. Rev. D53 (1996) 2958; J.M. Cline, K. Kainulainen 
and A.P Vischer, Phys. Rev. D54 (1996) 2451 

[10] M. Joyce, T. Prokopec and N. Turok, Phys. Rev. D53 (1996) 2930 

[11] P. Huet and A.E. Nelson, Phys. Rev. D53 (1996) 4578 

[12] P. Hernandez and N. Rius, Nucl. Phys. B495 (1997) 57 

[13] M. Carena, M. Quiros, A. Riotto, I. Vilja and C.E.M. Wagner, Nucl. Phys. B503 (1997) 387 

[14] M. Aoki, N. Oshimo and A. Sugamoto, Progr. Theor. Phys. 98 (1997) 1179; ibid. 98 (1997) 1325 

[15] J.M. Cline, M. Joyce, K. Kainulainen, Phys. Lett. B417 (1998) 79 

[16] A. Riotto, Nucl. Phys. B518 (1998) 339; Phys. Rev. D58 (1998) 095009 



15 



K. Farakos, K. Kajantie, K. Rummukainen and M.E. Shaposhnikov, Phys. Lett. B336 (1994) 
494; Z. Fodor, J.Hein, K. Jansen, A. Jaster and I. Montvay, Nucl. Phys. B439 (1995) 147; K. 
Kajantie, M. Laine, K. Rummukainen and M.E. Shaposhnikov, Nucl. Phys. B466 (1996) 189 

M. Carena, M. Quiros and C.E.M. Wagner, Phys. Lett. B380 (1996) 81; Nucl. Phys. B524 (1998) 
3 

D. Delepine, J.M. Gerard, R. Gonzalez-Felipe and J. Weyers, Phys. Lett. B386 (1996) 183 
J.R. Espinosa, Nucl. Phys. B475 (1996) 273 
M. Laine, Nucl. Phys. B481 (1996) 43 

J.M. Cline and K. Kainulainen, Nucl. Phys. B482 (1996) 73; Nucl. Phys. B510 (1998) 88 

J.M. Moreno, D.H. Oaknin and M. Quiros, Nucl. Phys. B483 (1997) 267; Phys. Lett. B395 (1997) 
234 

M. Losada, Phys. Rev. D56 (1997) 2893; G. Farrar and M. Losada, Phys. Lett. B406 (1997) 60; 
M. Losada, |hep-ph/ 99054411 

D. Bodeker, P. John, M. Laine and M.G. Schmidt, Nucl. Phys. B497 (1997) 387 
B. de Carlos and J.R. Espinosa, Nucl. Phys. B503 (1997) 24 

M. Laine and K. Rummukainen, Phys. Rev. Lett. 80 (1998) 5259; Nucl. Phys. B535 (1998) 423 

J.M. Cline and G.D. Moore, Phys. Rev. Lett. 81 (1998) 3315 

J.M. Moreno, M. Quiros and M. Seco, Nucl. Phys. B526 (1998) 489 

K. Kainulainen, talk given at COSMO-99, International Workshop on Particle Physics and the 
Early Universe, September 27 - October 2, 1999, Trieste, Italy; J. Cline and K. Kainulainen, work 
in progress. 

P. Elmfors, K. Enqvist, A. Riotto and I. Vilja, Phys. Lett. B452 (1999) 279 

A.G. Cohen, D.B. Kaplan and A.E. Nelson, Phys. Lett. B336 (1994) 41; D. Comelli, M. Pietroni 
and A. Riotto, Astropart. Phys. 4 (1995) 71 ||hep-ph/940^3"6T| 

S. Yu. Khlebnikov, Phys. Lett. B300 (1993) 376; A.G. Cohen, D.B. Kaplan and A.E. Nelson, 
Phys. Lett. B294 (1992) 57; J.M. Cline and K. Kainulainen, Phys. Lett. B356 (1995) 19 

D. Bodeker, Phys. Lett. B426 (1998) 351 

G.D. Moore and N. Turok, Phys. Rev. D56 (1997) 6533 



16 



